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Abstract 

In earlier work we provided the first evidence that the collapse, or coil-globule, transi- 
tion of an isolated polymer in solution can be seen in a four-dimensional model. Here we 
investigate, via Monte Carlo simulations, the canonical lattice model of polymer collapse, 
namely interacting self-avoiding walks, to show that it not only has a distinct collapse tran- 
sition at finite temperature but that for any finite polymer length this collapse has many 
characteristics of a rounded first-order phase transition. However, we also show that there 
exists a '0-point' where the polymer behaves in a simple Gaussian manner (which is a critical 
state), to which these finite-size transition temperatures approach as the polymer length is 
increased. The resolution of these seemingly incompatible conclusions involves the argument 
that the first-order-like rounded transition is scaled away in the thermodynamic limit to 
leave a mean-field second-order transition. Essentially this happens because the finite-size 
shift of the transition is asymptotically much larger than the width of the pseudo-transition 
and the latent heat decays to zero (algebraically) with polymer length. This scenario can be 
inferred from the application of the theory of Lifshitz, Grosberg and Khokhlov (based upon 
the framework of Lifshitz) to four dimensions: the conclusions of which were written down 
some time ago by Khokhlov. In fact it is precisely above the upper critical dimension, which 
is 3 for this problem, that the theory of Lifshitz may be quantitatively applicable to polymer 
collapse. 
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1 Introduction 



The collapse or coil-globule transition of an isolated polymer in solution has been studied by a 
variety of different theoretical approaches over the past 50 years ranging from phenomenological 
arguments, field theoretic renormalisation group approaches, continuum path integrals and the 
analysis of discrete lattice walks (l[ . Application and testing of these theories has mainly been 
confined to the "physical" dimensions of two and three. However, the phase transition of polymer 
collapse has been long believed to have an upper critical dimension of three {d u = 3) and 
so the differences between the predictions of many of the theories in that dimension lie in 
subtle logarithmic factors that are difficult to ascertain numerically In two dimensions the 
field theoretic (excluding conformal field theories) and continuum models have not given exact 
answers and cannot be compared to the conjectured || and numerically resolved Q values of 
universal quantities. Until recently || four dimensions has largely been ignored but here we 
argue that not only may four-dimensional studies be important to delineate which theoretical 
descriptions are valid, but that the collapse transition in four dimensions has some intriguing 
features of general interest in the field of phase transition and critical phenomena in statistical 
mechanics. 

An isolated polymer in solution is usually considered to be in one of three states depending on 
the strength of the inter-monomer interactions which are mediated by the solvent molecules and 
can be controlled via the temperature T. At high temperatures and in so called "good solvents" a 
polymer chain is expected to be in a swollen phase (swollen coil) relative to a reference Gaussian 
state so that the average size of the polymer scales with chain length algebraically faster than 
it would if it were behaving as a random walk. At low temperatures or in poor solvents the 
polymer is expected to be in a collapsed globular form with a macroscopic density inside the 
polymer. This implies an average size that scales slower than a random walk. Between these 
two states there is expected to be a second-order phase transition (sharp in the infinite chain 
length limit). 

The standard description of the collapse transition is a tricritical point related to the n — > 
limit of the (^ 2 ) 2 -(0 2 ) 3 O(n) field theory || |7], §]. One might then expect that above the upper 
critical dimension {d u = 3) some type of self-consistent mean-field theory based upon a suitable 
tricritical Landau-Ginzberg Hamiltonian [|| would give a full description of the transition, and 
hence conclude that in all dimensions d > 3 there is a collapse transition from a swollen state to 
the globular state with classical tricritical behaviour. There have been various other mean-field 
type approaches to this problem though their conclusions in three dimensions are similar || |l(]] 
to each other. 

The application of the mean-field theory of a tricritical point to polymer collapse predicts 
that at the transition point the polymer behaves as a random walk [y = 1/2), and this point has 
been known as the #-point (the #-point was originally defined as the point where the second virial 
coefficient of a dilute solution of polymers is zero, though it is expected that these definitions 
are equivalent). Thermo dynamically, one expects a weak transition with a jump in the specific 



heat a = (note that the thermodynamic polymer exponent a is related to the shift exponent 
ip = 2 — a in tricritical theory ||, itself not to be confused with the polymer theory finite-size 
scaling shift exponent). For finite polymer length N there is no sharp transition for an isolated 
polymer (unless one examines a macroscopic number of such polymers) and so this mean-field 
transition is rounded. In three dimensions the application of various self-consistent mean-field 
like approaches leads to the prediction that the second-order transition is rounded and shifted on 
the same scale of N~ 1 / 2 , that is, the crossover exponent <f> is 1/2, though strictly the power laws 
involved are modified via renormalisation group arguments [11, |l!| by confluent logarithms. (In 
particular note that it is predicted that in three dimensions the specific heat should be divergent 
logarithmically.) In four and higher dimensions no confluent logarithms should be present and 
one may expect pure mean-field behaviour with a crossover exponent of 1/2 (<pt is the relevant 
tricritical exponent here ||). 

On the other hand, Domb [|13| suggested some time ago that polymer collapse may be a first- 
order transition in three dimensions and the analysis by de Gennes || of the three-dimensional 
case of a suitably extended (Flory-type) self-consistent mean-field approach predicted in some 
parameter regions a first-order transition: this was superseded by his renormalisation group 
approach H In contrast some time ago there was the conjecture that the collapse transition 
disappears altogether above three dimensions, at least at finite temperature [El|. For d > 3 Sokal 



[16] has also pointed out that the field theoretic/Edwards model approaches have difficulties: in 
fact, if one analyses the Edwards model one finds the crossover exponent is given by <p = 2 — d/2, 
which for d = 4 gives = 0! In passing we note here that the same analysis predicts the shift of 
the #-point, defined say via the universal ratio of the radius of gyration to the end-to-end distance 
equalling its Gaussian value, should scale as N~^ 2 ~ 1 ' which has (polymer) shift exponent 1 
in dimension d = 4. This difference between the shift and the crossover exponent implies that 
strict crossover scaling has broken down. Of course, the theoretical fact that the swollen phase 
should also be Gaussian for d > 4 does raise the suspicion that the analysis of the Edwards 
model for polymer collapse may be subtle for d > 3. 

As a first attempt to explore the issues raised above we recently || considered the problem 
of interacting self-avoiding trails on the four-dimensional hyper-cubic lattice with a special set 
of Boltzmann weights as generated by a kinetic growth algorithm. Interacting self-avoiding 
trails are a candidate lattice model for polymer collapse. They are defined to be lattice paths 
such that each bond of the lattice may either be unoccupied or occupied by a single bond of 
the path, though sites can be multiply occupied. Attractive interactions are associated with 
those multiply visited sites so that the strength of this interaction drives a collapse. The study 
had several virtues. Firstly, since trails are allowed to intersect, but still possess excluded 
volume, we were able to grow configurations without the normal attrition that appears in growing 
self-avoiding walks through the "trapping" of the growth in dense sections of the polymer. 
Secondly, in lower dimensions it had been seen that whenever a kinetic growth algorithm was 
un-hindered by trapping it mapped precisely onto the collapse transition point of the model 



[17, |l^, In H we indeed found a set of Boltzmann weights where the model appears 

to have the Gaussian characteristics of a #-point. The two drawbacks of this approach are, 
firstly, that the simulations cannot effectively be extended away from the special temperature 
and, secondly, that the canonical lattice model of polymer collapse is rather self-avoiding walks 
(SAW) interacting via nearest-neighbour site (monomer) attraction. This canonical model is 
known as interacting self-avoiding walks, hence ISAW. A new algorithm, known as PERM, 
for the Monte Carlo simulation of ISAW (among other things) has been recently developed by 
Grassberger and collaborators [20, [H], |2^[. This is essentially a kinetic growth strategy that adds 
clever enhancements to simultaneously allow a wide range of temperatures to be accessed and for 
the attrition of samples through trapping to be lessened (see section [| for a fuller description). 
In this paper we have simulated ISAW on the four-dimensional hyper-cubic lattice using the 
PERM algorithm over a wide range of temperatures. 

Now, to begin, our results suggest that there is indeed a collapse transition in four dimensions 
at a finite temperature. However, the character of that transition is particularly intriguing! We 
find a distinct double peak distribution for the internal energy far below a point which we clearly 
identify as a candidate #-point. This double peak distribution becomes more pronounced as the 
chain length is increased. This would seem to suggest a first-order transition. If this was the 
case there would be a delta function peak forming in the specific heat but we find that while a 
peak is indeed forming it seems not to be growing linearly with the size of the polymer. More 
importantly, the location of a distinct 0-point is incompatible with a first-order transition if 
there is only one collapse! However there is a theoretical framework (whose conclusions are 
suitably extended here) that is consistent with the evidence we present. This framework was 
explained in a paper by Khokhlov [23] who applied the mean-field approach of Lifshitz, Grosberg 
and Khokhlov (LGK) |To| , 24, 25 to arbitrary dimensions^. The LGK theory is based on a 
phenomenological free energy in which the competition between a bulk free energy of a dense 



globule and its surface tension drive the transition. Until recently [26] the consequences of this 



surface free energy were largely ignored in the polymer literature. Its effect on the scaling form 



of the finite-size partition function was argued and confirmed [26, 27, pq] . We shall refer to 
the LGK theory as applied to dimensions four and above as KLG to distinguish it from the 
original three-dimensional work of Lifshitz, Grosberg and Khokhlov [p4| , |25| . 

Hence the major conclusions of our work are that the finite-size character of the coil- globule 
transition in four dimensions is first-order despite the thermodynamic limit being probably 
adequately described by mean-field tricritical behaviour. The only alternative conclusion from 
our data is that the transition is truly first-order and our finding of a 0-state is fortuitous. The 
whole theory of crossover scaling for this transition needs to be reworked. This curious state 
of affairs where a second-order transition looks distinctly first-order may be of interest in other 
physical situations where mean-field theory is used to describe thermodynamics. 



x We warn the reader that the abstract (and parts of the conclusions) of the paper [£3| may be misleading as 
it reads that "for d > 3 the coil-globule transition is of first-order". Moreover the application of this mean- field 
theory to d < 3 is inappropriate and its conclusions have now been superseded. 



The layout of the paper is as follows. In the next section we define the model we consider and 
review the generally expected behaviour of the quantities we have calculated in four dimensions. 
Then in section 3 we explain the results of KLG theory as applied to four dimensions. In section 
4 we explain our Monte Carlo approach, PERM, and finally in section 5 we carefully describe 
the numerical results of our simulations and how well they conform to the theory of KLG. 



2 The ISAW model and a review of basic scaling results 

The interacting self-avoiding walk model is the canonical lattice model of the coil-globule tran- 
sition and has been long studied in two and three dimensions. Here we shall consider the 
four-dimensional hyper-cubic lattice (coordination number 8). The monomers are imagined to 
be sitting on the sites of the lattice and a self-avoiding path of such sites form the polymer. The 
self-avoidance means that no two monomers can sit at the same site of the lattice. 

The partition function of the self-interacting self-avoiding walk model (ISAW) is given by 

Z N {u) = ]T , (2.1) 

where the sum is over the set of all self-avoiding walks Q.^ of length N steps (N + 1 monomers) 
with one end at some fixed origin and m(<p) is the number of non-consecutive nearest-neighbour 
monomers for a given walk ip. The Boltzmann weight ui = e^ e is associated with a nearest- 
neighbour contact of energy — e so that to > 1.0 for attractive interactions. We define a reduced 
finite-size free energy per step kn(w) as 

K N {uj) = ^\ogZ N {u). (2.2) 

The usual free energy is related to this by —f3F^ = Nkn(uj). 

The average of any quantity Q over the ensemble set of allowed paths £In of length N is 
given generically by 

WhM _ E^W"*" . (2 . 3) 

We define a normalised finite-size internal energy per step by 

U N (u) = ^ , (2.4) 
and a normalised finite-size specific heat per step by 

These quantities are related in the usual way to the reduced free energy via Un = 8kn /dlogu 
and Cn = 9Un /dlogu. 

The thermodynamic limit in this problem is given by the limit ./V — > oo so that the thermo- 
dynamic free energy per step foo{uj) is given by 

-/5/ooM = «ooM = lim k n {u) . (2.6) 



This quantity determines the partition function asymptotics, i.e. Z^(u) grows to leading order 
exponentially as ^(oj) n with /j,(uj) = e Kac ^. 

In our simulations we calculated two measures of the polymer's average size. Firstly, specify- 
ing a walk by the sequence of position vectors ro, ri, r^r the average mean-square end-to-end 
distance is 

{R 2 e ) N = ((rjy - r ) • (rjv - r )) • (2.7) 
We shall use the symbol R 2 n to be equivalent to 

Rl N {u) = (R 2 e ) N . (2.8) 

The mean-square distance of a monomer from the endpoint ro is given by 

1 N 

(On = j^T[ - r o) • ( r * - r o)> • (2.9) 



8=0 



Again we define 



We also define the ratio 



R 2 m M = (R 2 m ) N . (2.10) 



B N (u) = |K (2.11) 

e,N 

which should have a universal limit in each critical phase of the model. 

Now let us assume for a moment that there is a single collapse transition at some value of 
lo and let us explore the (four-dimensional) behaviour we might expect from the above defined 
quantities in each of the phases. As discussed above, the basic physics of the coil-globule 
(collapse) transition can be understood by the consideration of the average size of the polymer, 
-R/v, either R e< N or R m ,N, as a function of length N in each of the phases, so let us consider this 
first. Generally one always expects that 

R 2 N ~ a(oj) N 2v as N ^ oo (2.12) 

for any fixed temperature. In four dimensions at infinite temperature, uj = 1, it has been 
predicted [^9| that 

R 2 N ~a + N (log(AO) 1 / 4 . (2.13) 

If there does exist a collapse transition then one would expect that this scaling extends (with a 
constant a + that depends on temperature) down to the transition point. In the collapsed phase 
the polymer is expected to assume a dense configuration on average and hence the globular value 
of the radius-of-gyration exponent is v g = 1/d = 1/4 || with 

R 2 N ~ aT{u) N 1/2 . (2.14) 

Finally at some finite transition temperature 1.0 < Ut < oo a Gaussian scaling of the radius of 
gyration should occur, that is 

R 2 N ~ a N, (2.15) 



so that u t = 1/2. This Gaussian scaling is often used (theoretically at least) to define the 0-point 
u) = ujg of an isolated polymer so that uj t = ujq. The universal ratio is expected to converge 
to the value -Boo = 1/2 both in the swollen phase and at ujg. However, one would expect slow 
logarithmic corrections for oj < ujg and algebraic corrections at loq. For u) > uig the phase is 
no longer expected to be critical and so is no longer universal and may be a non-constant 
function of u. 

One can also consider the scaling of the partition function in each of the regimes, given that 
there is a transition. For high temperatures 1.0 < u < u$ one expects the infinite temperature 



behaviour, which is [29| 



Z N ~ b + (w) n(u) N (log iV) 1/4 , (2.16) 
while at low temperatures [26] one expects asymptotics of the form 

Z N ~b-(u) fi{uj) N v s {oj) Ni/i N 9 (2.17) 

where [i s is related to the surface free energy of the polymer globule and the exponent g need 
not be universal (we only write it for completeness of the asymptotic form). For uj = ujq one 
expects 

Z N ~ b e h(lu 9 ) N (2.18) 

as a reflection of Gaussian behaviour. 

In the thermodynamic limit the thermodynamic functions foo(^), t^oo(w) and Coo(u;) are all 
expected to be analytic functions of u except at ojq. By using the correspondence to the tricritical 
model Q the mean field theory would imply that the specific heat had a jump discontinuity 
at ujq since a = 2 — ip = 0. Of course, for finite N there is no sharp transition for an isolated 
polymer (unless one examines a macroscopic number of such polymers). Around the collapse 
temperature in three dimensions the finite-size corrections are expected to take on a cross-over 
scaling form || so that 

R% ~ a 9 N1Z{{T -Te)N^) (2.19) 

with = 1/2 but this form implies that the transition is rounded and shifted on the same scale 
of and so its applicability to four dimensions needs careful thought. 



3 The theory of Khokhlov, Lifshitz, and Grosberg (KLG) 

We now review the results of Khokhlov [^] paying special attention to the predicted behaviour 
of the quantities we have calculated. The theory originally proposed by Lifshitz for the general 
mean-field description of the globular state, extended by Lifshitz, Grosberg and Khokhlov to 
fully describe the transition, and finally applied to four dimensions by Khokhlov |23|] is based 
on several phenomenological mean- field assumptions. Firstly, there exists a state where the 
excluded volume property of long chain molecules is exactly cancelled by the attractive interac- 
tions between parts of the polymer as mediated by the solvent. This is the #-state. Secondly, 



when the attraction becomes even stronger there eventuates a globular state where the polymer 
behaves as a liquid drop. The results of the theory are based on a phenomenological free energy 
of that globular state relative to the free energy of the pure Gaussian state of the #-point at 
Tg. Hence the condition applied to find the finite-size position of the transition is to equate the 
relative free energy to zero. The relative free energy is given as 

Fn = Fbulk + Fsurface (3.20) 

where the F^ u ik and F sur f ace are given in terms of the second and third virial coefficients, the 
length of the chains, and the linear size of the polymer found from the globular density. In 
particular both the bulk and surface free energies are proportional to the square of the second 
virial coefficient. It is assumed that on approaching the #-point the second virial coefficient goes 
to zero linearly with temperature while the third virial coefficient remains non-zero. Note that 
this implies a quadratic dependence of the bulk free energy on the distance to the 0-point. Since 
the free energy has exponent 2 — a this implies an exponent a = (assuming that this part of the 
free energy is singular). Therefore a second-order phase transition occurs in the thermodynamic 
limit. 

It is further assumed that the density in the globule is proportional to the second virial 
coefficient and hence also goes to zero linearly with temperature on approaching the #-point (/? = 



1). By applying the condition Fn = 0, Khokhlov [ 23 1 finds a finite-size transition temperature 
T Ci jv that approaches the ^-temperature from below as iV -1 / 3 . Applying this result to the ISAW 
model gives 

s 

u c ,N -uj e ~ (3.21) 

for some constant s. That is, the polymer collapse shift exponent is 1/3. 

Now Khokhlov [23] found that the free energy (for T < Tg) can be rewritten in terms of T c< n 



as 



-NT (T e — T) 1 



(- 


T$ — T c> n 


3/4\ 




T e -T 





(3.22) 



Actually KLG's unwritten assumption is that Fn is the minimum of zero and the right-hand 



side of ( |3.22[) so it is zero for Tg < T < T Cj n- Khokhlov 2J| deduces that the width of the 
transition AT at finite N can be found from this and scales as iV~ 2 / 3 . Hence the ISAW model 
should have a transition width Aw that scales as 

w 

Aw ~^ ( 3 - 23 ) 

for some constant w. That is, the polymer collapse crossover exponent is 2/3. Hence note that 
the size of the crossover region is asymptotically small relative to the shift of the transition, and 



that crossover forms such as (2.19) may not be useful. 

Even though there can be no sharp transition for a single polymer of finite length N the 
theory can describe the nature of the rounded transition by considering the difference between 
the density of the globular state p g at T c> n and that of the coil state p c at the same temperature. 



This is 

P g (T c ,N) ~ Pc(T c ,n) 



Pc{T t 



c.N) 



N 2 / 3 , (3.24) 



which diverges as N becomes large. Hence Khokhlov [23 concluded that 'the coil-globule tran- 
sition is first-order', though we now interpret this to mean that the finite-size corrections to 
the thermodynamic second-order transition are first-order likef]. However, both p g {T c ^) and 
Pc{T c ,n) tend to zero as N —* oo and it is simply that p g (T c ^) tends to zero asymptotically 
slower than p c (T c ^) that makes the relative difference quoted above diverge. Noting that 
p = N /R A and R c ~ iV 1 / 2 , the above equation can be used to deduce the scaling of Rn at T Cj tv 
as 

Rn(T CjN ) ~ a c N^ 3 . (3.25) 

Hence we define an effective radius-of-gyration exponent v c = 1/3 for the scaling of the size 
of the polymer when following the finite-size transition temperatures. Note that this exponent 
value obeys vq = 1/2 > u c > v g = 1/4. 

Following the the work of Lifshitz, Grosberg and Khokhlov [^] one can also calculate the 
change in the internal energy over the crossover width of the transition AT as the latent heat 
(or "heat of the transition") by using expression (3.22): 



The corresponding height of the peak in the specific heat is 

C N (T C>N ) ~ h c iV 1 / 3 . (3.27) 

So to summarise the picture so far, the theory predicts a thermodynamic second-order tran- 
sition at a #-point with a jump in the specific heat. For finite polymer length this transition 
is shifted below the #-point by a temperature of the order of 0(N^ 1 ^ 3 ) with the width of the 
transition of the order of 0(iV~ 2 / 3 ). Over this width there is a rapid change in the internal 
energy that scales as 0(A r ~ 1 / 3 ): the important point here of course is that this tends to zero 
for infinite length so the effect of the peak in the specific heat is scaled away for iV large, leav- 
ing a finite jump in the thermodynamic limit. To understand this further let us consider the 
distribution of internal energy as a function of temperature and length. For any temperature 
above the 0-point and well below T c ^ one expects the distribution of internal energy to look 
like a single peaked distribution centred close to the thermodynamic limit value: a Gaussian 
distribution is expected around the peak with variance 0(N~ 1 ^ 2 ). In fact, this picture should 
be valid for all temperatures outside the range [T C) jv — 0(N^ 2 ^),T c ^ + 0(iV -2 / 3 )]. When we 
enter this region we expect to see a double peaked distribution as in a first-order transition 
region. For any temperature in this region there should be two peaks in the internal energy 
distribution separated by a gap SU of the order of 5U ~ AU oc 0(N^ 1 / 3 ). Each peak is of 

2 We point the terminology of Khokhlov was presumably that explained in Section I.C.2 of [[^5| but may be 
misleading to the modern reader. 



Gaussian type with individual variances again of the order of 0(N~ 1 / 2 ). Hence as N increases 
the peaks will become more and more distinct and relatively sharper but the peak positions 
will be getting closer together. We refer to this scenario as a pseudo-first-order transition or, 
more correctly, as first-order- like finite-size corrections to a second-order phase transition. If 
there were a real first-order transition then the distance between the peaks should converge to 
a non-zero constant. 



4 Perm 

We have simulated IS AW using the Pruned-Enriched Rosenbluth Method (PERM), a recently 
proposed generalisation of a simple kinetic growth algorithm PERM builds upon the 

Rosenbluth-Rosenbluth method [3C], in which walks are generated by simply growing an existing 
walk kinetically, i.e. by choosing the next step with equal probability from all possible accessible 
lattice sites. Eventually, a walk generated thus gets trapped in a configuration in which it cannot 
be continued, leading generally to an exponential "attrition". Moreover, in order to simulate 
ISAW at a particular temperature one needs to re-weight the kinetically grown samples in such 
a way that the generated sample is usually dominated by a few configurations which carry large 
weight after the re-weighting. 

In order to overcome these obstacles, PERM uses a combination of enrichment and pruning 
strategies to generate walks whose weights are largely distributed around the expected peak of 
the distribution. On the one hand, if the weight of a configuration becomes too small, the con- 
figuration gets pruned probabilistically and the weight adjusted correspondingly. Alternatively, 
if the weight of a configuration becomes too large, copies of the walk are made and the respective 
weights reduced accordingly. While this does not eliminate trapping, it is generally sufficient 
to compensate for it: trapping occurs when the end of the walk is in an area of high density, 
which in turn increases the likelihood of enrichment. The algorithm can be implemented in a 
self-tuning way by choosing dynamically adjusted upper and lower threshold values to control 
pruning and enrichment rates. 

It is plausible that the algorithm works best at temperatures in which the thermal distri- 
bution is close to the distribution of walks generated by kinetic growth. In sufficiently large 
dimensions, this temperature should be quite close to the ^-temperature, so that the algorithm 
is expected to be well suited to the study of polymer collapse. As mentioned in the introduc- 
tion, there are lattice models of interacting polymers for which there is an exact mapping of 
the corresponding kinetic growth models to their respective /9-points, in which case a PERM 
simulation at the ^-temperature reduces to simple kinetic growth. In fact, it turns out that the 
algorithm performs well over a whole range of temperatures covering all of the swollen phase 
and the scaling region around the collapse transition. However, we find that the performance of 
PERM in the collapsed phase is far less satisfactory. 

The guiding principle for any choice of implementation should be the observation [EG] that 



the algorithm essentially produces a random walk in chain length with reflecting boundaries 
at and N max . Considered in such a way, the algorithm performs best if this random walk is 
unbiased and if the associated diffusion coefficient is large. To eliminate bias, pruning, trapping 
and enrichment rates have to compensate each other. To maximise the diffusion coefficient, the 
pruning and trapping rates have to be as small as possible. The choice of pruning and enrichment 
thresholds needs to take both into account. (In contrast with the original work on PERM, where 
trapping was viewed as a special case of pruning, we find it instructive to distinguish between 
these two effects: trapping is unavoidable due to the geometry of the lattice, whereas pruning 
is done optionally to adjust weights.) 

In our implementation, we chose upper and lower thresholds W u and W l proportional to the 
current estimate of the average weight of a walk at length N, (Zn)/sn, where sn is the number 
of generated samples at length N, and (Zn) is the current estimate of the partition function at 
length N. That is to say, 

W% = c u N (Z N )/s N , W l N = c l N (Z N )/s N . (4.28) 

In order to enforce an even sample size distribution we allow for dynamic adjustment of cfj and 
c l N . Thus, if for example at any particular length we have an excess of pruning, the algorithm 
"relents" and increases both Cjy and c l N in order to reduce pruning and enhance enrichment, 
keeping the quotient of the thresholds Q = c^/c l N constant. To stabilise the dynamic adjust- 
ment, we enforce Cjy > c^ in and c l N < c l max . After some initial experimentation, we chose 
c min = ^ and c l max = 1/2, which leaves us with the threshold quotient Q as the sole adjustable 
parameter. (We also experimented with a dynamic length-dependent Qn, but the dynamic ad- 
justment seemed to be too unstable to pursue this avenue further.) For each run, we attempted 
to choose the smallest threshold quotient Q for which we could obtain an even sample size 
distribution. 

The disadvantage of PERM is that due to the enrichment the generated data is not inde- 
pendent. All the data generated during one "tour", i.e. between two successive returns of the 
algorithm to length 0, is correlated. Therefore, we keep track of the statistics of tour sizes t 
to get a rough idea of the quality of the data. In our statistical evaluation we use (somewhat 
arbitrarily) the quotient of sat and \J (t 2 ) as a measure of an effective independent sample size. 
This is correct as long as the tour sizes don't fluctuate too strongly, and, more importantly, 
as long as individual tours explore the sample space evenly. When simulating in the collapsed 
phase, both of these assumptions break down, and the sample is dominated by few huge tours. 
Moreover, the pruning and enrichment rates become so large that the efficiency of the algorithm 
is significantly decreased. 

For further details of the algorithm and suggestions of various other improvements such as 



Markovian anticipation we refer to [21, 22 1. 



5 Results 



We simulated ISAW on a 4-dimensional hyper-cubic lattice using PERM with N max set to 1024, 
2048, 4096, 8192, and 16384 at values of uj ranging from 1.0 to 1.4 {N max = 1024) to 1.2175 to 
1.2225 (N max = 16384) and the threshold quotient Q ranging from 10 to 160. At each fixed uj, 
we generated 10 7 walks. To illustrate the computational effort, the generation of a sample of 
size 10 7 at length N max = 16384 took about 2 weeks CPU time on a 600 MHz DEC Alpha. 

We computed statistics for N and R^n^ the partition function Zn, the internal energy 
Un and specific heat Cat. Moreover, we generated the distribution of the number of interactions 
at N max . The distributions obtained at various temperatures were then combined using the 



multiple histogram method 31]. 

In discussing our findings on the nature of this polymer collapse transition, it is natural to 
first present the change of size of the polymer as the interaction strength changes. In Figure [l], 
we display the mean-squared end-to-end distance normalised by walk length, R^, N /N for lengths 
1024 up to 16384 (we have analogous data for the quantity R 2 m N /N). As discussed, one expects 
this quantity to increase logarithmically in A" in the swollen regime, approach a constant at a 
random-walk like 0-point and decrease in iV as a power law in the collapsed regime. As can 
be seen clearly from Figure |l|, there is indeed a transition from a swollen region, where this 
parameter increases with N to a collapsed region, where the value has dropped sharply. In the 
transition region, however, two phenomena can be noticed. On the one hand, around uj = 1.18 
the quantity R^ N /N approaches a constant, which is indicative of #-point behaviour. On the 
other hand, the collapse occurs in a region which is well separated from this 0-point. With 
increasing polymer length, the region where the collapse occurs approaches the #-point, but 
simultaneously sharpens so strongly that it remains well separated from it. 

In the swollen phase, our results are in correspondence with the logarithmic corrections seen 
by Grassberger see Figure ||. As in that paper, we observe that R^ N grows as N(logN) c , 
albeit not with the exponent predicted by field theory. At uj = 1.0 we find c = 0.30. This value 
shifts to c = 0.22 at uj = 1.1, indicating the presence of strong temperature-dependent correction 
terms. 

Near the suspected #-point, we extended our simulations to walks of length 32768. Figure 
^ shows a plot of R^ N /N versus 1/N for values of uj between 1.180 and 1.184. At uj = 1.182 
we have indeed an asymptotically linear dependence of R^n on N. Moreover, at uj = 1.182 
we estimate from our data = RmN/^-eN = 0.500(1), which is also indicative of Gaussian 
behaviour. 

As was already seen from Figure |], the collapse happens very rapidly as uj increases. An 
alternative way of visualising this is to consider how the size of the polymer changes at fixed uj in 
the collapsed phase as the length A" increases. As shown in Figure || for uj = 1.4, R? e N changes 
non-monotonically in AH After an initial increase, the size of the polymer actually shrinks 
around A" = 250 as it undergoes collapse^] corresponding to a rapid increase of the density. For 
3 Note that this lack of monotonicity is not an indication of a first-order transition. A similar non-monotonous 



large enough N, we expect to see the true collapsed behaviour, i.e. R 2 N growing again as N 1 / 2 , 
but this regime is beyond the reach of our PERM simulations on current computer hardware. 

The swollen phase and the 0-point behaviour can also be clearly identified from the free- 
energy scaling. In the swollen phase we find again the same behaviour as [32]. Figure || shows 
Zj\f//j N for uj = 1.0 and uj = 1.1. For each uj, this quantity is plotted with three values of /i 
which differ in the sixth digit, showing both the accuracy in the estimation of the free energy 
and the presence of logarithmic corrections. At uj = 1, we estimate fj>(l) = usaw = 6.77404(2). 
This can be compared to earlier estimates of 6.7720(5) 6.774(5) @, and 6.77404(3) |2| 

In the 0-region, a similar analysis shows that here Zjy scales as /j, n with weak 1/N corrections. 
Figure |6| shows Zn/h n plotted versus 1/N, with respective values of [i obtained in a similar 
fashion to the one shown in Figure ||. We estimate the #-point to be ujq = 1.182(1) and /xg = 
7.011(2). (At fixed uj, the accuracy is of course higher: for ui = 1.182, we estimate fj, = 
7.0117675(5).) 

In the collapsed region, one expects the finite-size free energy to have a strong correction 
term of the order N~ l / A due to surface effects. Figure ^ shows this for ui = 1.4. As argued 
above, the globule collapses when the length is above N = 250, and we notice here the onset 
of a corresponding strong change in the behaviour of the finite-size free energy around this 
length (A^- 1 / 4 ~ 0.25). Even though we cannot simulate long enough chain lengths to clearly 
determine the precise nature of the correction term, our data is certainly compatible with a 
JV-V4 correction for N~ 1 / 4 < 0.2. 

In order to study the collapse transition more closely, we now focus our attention on the 
internal energy and specific heat. As can be seen from Figures || and ^, the internal energy 
increases rapidly over a small temperature interval with a corresponding diverging specific heat. 
As N increases, the transition becomes sharper and stays well separated from the #-point, even 
though the location of the transition approaches the #-point slowly. 

The scaling of the shift of the transition towards the #-point uj c n — ujg and the sharpening 



of the transition width Aw are both shown in Figure 1C . Here, we defined the location of the 
collapse transition by the location of the specific heat peak, and the width of the transition is 
given by the interval in which the specific heat is greater or equal to half the value of the peak 
height. Expecting from the KLG theory that oj C) n — ujq scales as A r_1//3 and that Aw scales as 
A^ 2 / 3 , we plot both A rl / 3 (o; c jv — ujq) and AujN 2 ^ 3 versus N~ 2 / 3 which was chosen empirically. 
Both quantities can be seen to be asymptotic to constants: on the graph extrapolations give 



non-zero intercepts. Hence, Figure 1C shows that the KLG predictions are compatible with 
our simulations. We do note that the corrections to scaling for Auj are much larger than for 

U C ,N — ^0. 

The character of the transition becomes apparent if one plots the internal energy density 
distribution (rescaled density of interactions) at the finite-size collapse transition, uj c ^. Figure 



11 shows the emergence of a bimodal distribution. At length 2048 one sees a slight non-convexity, 



behaviour can be observed for ISAW in three dimensions, where the collapse transition is second-order 



which at length 16384 has evolved into a distribution dominated by two sharp and well-separated 
peaks. The values of the minima and maxima of the distribution are different by two orders of 
magnitude. 

It is instructive to study the transition by how this distribution changes over a large range 



of uj. Figure 12 shows this for N = 4096. One sees that there is not much of a change in 
the shape and location of the distribution between the non- interacting case, uj = 1, and the 
#-point, uj = 1.182. However, in a very small interval around the collapse transition, the density 
distribution changes dramatically as uj increases. The density distribution switches from a peak 
located around 0.3 to a peak located around 0.55, corresponding to a sudden change in the 
internal energy. In the collapsed phase, the width of the peak is much wider than in the swollen 
phase, implying a larger specific heat. It is this difference between the swollen and collapsed 
phases' specific heats that will eventually become the thermodynamic second order jump. 

The rapid first-order like switch between two peaks in the distribution becomes more pro- 
nounced at larger polymer lengths. At N = 16384, this 'switching' is shown in Figure |l^: over a 
range of uj of the order of 10~ 3 , a peak near uj = 0.3 disappears while another peak near uj = 0.5 
emerges. 

Returning to the scaling predictions from KLG theory, a suitably defined finite-size latent 
heat, AQ, should tend to zero as iV -1 / 3 in the thermodynamic limit. One possible measure 
of this latent heat is given by the product of specific heat peak Cn(^c,n) and specific heat 
width Aw, and another is given by the distance 5U of the peaks in the bimodal internal energy 



distribution. Figure |14| shows the behaviour of both of these quantities. One notices two things 
from this figure. Firstly, we are unable to confirm or deny the predicted scaling behaviour for 
AU and, secondly, even at length N = 2048 (iV -1 / 3 ~ 0.08) there is considerable discrepancy 
between the two quantities, so that it is not surprising that one cannot discern a clear scaling 
behaviour. The explanation for the discrepancy between the two quantities as well as of the 
difficulty of observing the predicted scaling behaviour is of course that in order to observe the 
asymptotic behaviour the two peaks in the histogram have to be well separated and distinct, 
and we see from Figure |ll] that this is only the case when N is of the order of 10 4 . This explains 
why we are unable to find a value for the exponent related to the divergence of the specific heat 



consistent with the rest of our theoretical picture. We do concede that Figure 14 alone could be 
used to argue for the existence of a real first-order transition in the thermodynamic limit, but 
we believe the rest of our data and other theoretical facts provide a more consistent picture. 

In conclusion, our ISAW simulations elucidate the structure of the polymer collapse tran- 
sition in four dimensions. We show conclusively that there is indeed a collapse transition at 
a finite temperature. Secondly, we find evidence for a ^-temperature at which the polymer 
is well approximated by Gaussian behaviour as well as for a collapse transition which is well 
separated from the #-point. The collapse transition shows many first-order like features, such 
as a bimodal distribution in the internal energy. An analysis of the scaling behaviour of this 



transition in the context of the theory of Lifshitz, Grosberg and Khokhlov p5|, 23] shows that 



a consistent interpretation of these findings is that of first-order like finite-size corrections to a 
thermodynamic second-order transition. We note that these findings are reminiscent of results 
for interacting self-avoiding trails on the diamond lattice, where a #-point was found in jl{J and 
subsequent simulations revealed a bimodal distribution in the internal energy density []35|]. In 
[19] it was concluded that the transition is second-order, whereas in |35|] the conclusion was that 



the transition is first-order. In the context of the findings presented here, it is tempting to expect 
a similar resolution of this apparent contradiction in terms of a pseudo-first-order transition. 
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Figure 1: R% N /N versus u for lengths N = 1024, 2048, 4096, 8192, and 16384. 
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Figure 2: R^ N /N versus N at uj = 1.0 and w = 1.1. The curves are fits to R^ N /N = alog(N + 
b) c over the shown range. 
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Figure 3: R\ N jN versus l/N in the 6-region: uj = 1.180, 1.181, 1.182, 1.183, 1.184 from top to 
bottom. 
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Figure 4: R 2 e N versus N for lo = 1.4. 
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Figure 5: Zjy / 'fi N versus N at to = 1.0 and u> = 1.1. At each value of u, three different values 
of are shown. We estimate //(1.0) = 6.77404(2), and /i(l.l) = 6.89699(2). The filled circles 
correspond to the central estimates of fi while open circles correspond to shifting fi by the error 
estimate quoted. 



1.22 



1.21 - 



1.20 - 



Z N /fi 



N 



1.19 



1.18 



1.17 





O 


o 


o 


o 


o 


o 




+ 


+ 


+ 


+ 


+ 




**** * 








* 


* 




Ax X 


X 


X 


X 


X 


X 


X - 












• 


• 








• 


• 






• 


• 


• 












1 




1 




1 1 





0.0000 0.0002 



0.0004 0.0006 
1/N 



0.0008 0.0010 



Figure 6: Z^/fi N versus 1/N using our best estimates of /i in the 6-region: u> = 1.180, 1.181, 
1.182, 1.183, 1.184 from top to bottom. 
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Figure 7: Finite-size free energy km versus N 1 / 4 for to = 1.4. 




Figure 9: Specific heat C N versus uj for lengths N = 1024, 2048, 4096, 8192, and 16384 from 
right to left respectively, using the multi-histogram method. 
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Figure 10: Scaling of the transition: shift and width of the collapse region. Shown are the scaling 
combinations 7V 1 / 3 (a; Ci iv — wq) an d N 2 / 3 Alu versus N~ 2 / 3 . 




Figure 11: Internal energy density distributions at uj C) n for 2048 and 16384. The more highly 
peaked distribution is associated with length 16384. 




Figure 13: Internal energy density distributions at uj = 1.2195 and 1.2210 for N = 16384. 
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Figure 14: Scaling of the latent heat AU: our two measures of AU, Cjv(o; C) jv)Aa; and peak 
distance 5U are plotted versus iV -1 / 3 . 



